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OPTIMIZATION  OP  SPACE  TBAJECTOBIBS 


S.  I.  Kirpichnikov. 

i 

the  works  of  this  section  are  dedicated  to  tke  development  of 
aetkods  of  the  aechaaics  of  optiaal  coatro liable  systems  as  applied 
to  prpbless  of  space  dynamics.  The  primary  number  of  works  is 
coaaeoted  with  the  dewelopoeat  of  aaalytlcal  methods  of  the  options 
select ioa  of  traasfer  orbits  io  a ceotral  field  aod  with  the 
approximate  solution  of  separate  model  problems^ 

the  work  of  7.  A.  Antonov  and  A.  S.  Shayrov,  placed  in  the 
section  of  brief  scientific  reports*  is  directly  adjacent  to  the 
articles  of  this  section. 


T 
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OPTIHJBATIOR  or  dsscbst  habbovbbs  op  satellite  bodole  pbob 

IBAB-PLAIBT  OBBIT 


Let  as  exanine  a space  vehicle#  noviag  aloqg  an  elliptical  orbit 
iq  the  gravitational  field  of  a spherically  syaaetric  planet.  The 
space  vehicle  consists  of  two  parts:  aodele  descendable  to  the  planet 
(nod she  i)  aad  orbital  nodale  (nodale  II)  <>  It  is  required  to 
constsact  an  optiaua#  with  respect  to  fuel  consumption,  descent 
aaneaver  of  sodale  I to  the  planet.  The  aaneaver  is  executed  with  the 
aid  of  a single  initial  pulse,  applied  to  nodule  I so  that  the 
orbital  nodule  continues  to  nove  along  the  initial  orbit,  and  the 
descent  nodule  transfers  to  the  approach  trajectory  to  the  planet. 

The  landing  naneuver  of  nodule  I is  terninated  after  entry  into  the 
dense  layers  of  the  ataosphere  and  their  passage.  The  planet  nay  not 
have  an  at nosphere,  then  the  naneuver  is  terninated  by  the  hard 
contact  of  nodule  1 with  its  surface.  In  the  latter  case  it  is 
necessary  everywhere  below  to  fornally  egaate  the  altitude  of  the 
at nosphere  to  sero. 

Let  us  agree  by  start  to  aean  the  aoaent  of  the  initial  pulse. 

By  finish  for  deterninacy  we  will  nean  the  landing  of  nodule  I on  the 
surface  of  the  planet.  However#  all  the  discussed  stops  being  valid 
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if  as  the  aoaeat  of  finish  we  select  sose  other  fixed  soaent  between 
entry  of  sodale  I into  dense  layers  of  the  ataosphere  and  its 
landing. 

The  initial  orbit  of  the  space  vehicle  and  the  trajectory  of  the 
descent  sodale  right  ap  to  its  entry  into  the  dense  layers  of  the 
atnosphere  are  considered  Keplerian,  where  there  are  considered  only 
elliptical  orbits  of  descent  with  straight  notion  with  respect  to  the 
initial  orbit.  It  is  assaaed  that  the  angular  range,  change  of 
altitude  and  the  tine  of  notion  on  the  segsent  of  passage  of  dense 
layers  of  the  atnosphere  to  the  finish  point  are  known  beforehand, 
and  that  this  segnent  is  located  in  planetocentric  plane,  passing 
through  the  velocity  vector  of  nodule  I during  its  entry  into  the 
dense  layers  of  the  ataosphere. 

Taking  into  account  the  Tsiolkovskiy  fcraula,  we  will  uiniaise 
the  characteristic  velocity  of  the  initial  pulse,  i.e. , the  aodulus 
of  pulse  change  of  velocity  of  nodule  I. 


In  this  work  the  foraulated  problea  is  igvestigated  with 
coaplicated  boundary  conditions.  lanely  limitations  are  introduced  on 
the  angle  of  entry  of  the  descent  nodule  and  on  the  distance  between 
the  apdules  at  the  aoaeat  of  finish.  Purtheraore,  there  is  considered 
the  cpadition  of  direct  visibility  between  the  nodules  and  nore 
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general  liaitation  on  the  co-altitude  of  aodule  n at  the  point  of 
finish,  and  there  can  be  taken  into  accoaat  the  liaitation  on  the 
velocity  of  aodale  I at  the  aoaent  of  eatery. 


following  the  aethod  proposed  by  Ting  Lu  £1]  it  is  easy  to  prove 
that  the  options  descent  orbits  of  the  aaaenvers  exaaiaed  here  in  all 
cases*  being  of  practical  interest,  should  lie  in  the  plane  of  the 
initial  orbit:  therefore  this  investigation  is  United  to  coplaner 
for natation  of  the  problea. 

Q1.  Batheaatical  F cumulation  of  the  problea  and  General 
Conclusions. 

In  the  plane  of  notion  let  os  in  trod  ace  polar  coordinates  r<  * 
with  the  origin  at  the  center  of  the  plaqet  so  that  the  direction  of 
positive  reading  of  angle  4 coincides  with  the  direction  of  notion 
along  the  initial  orbit. 

hloag  with  Kepler ian  eleaeats  of  orbits:  large  seniaxis  a, 

eccentricity  e and  angular  distance  e of  pericenter  froa  the  polar 

* 

axis  as  will  consider  eleaeats  p,  g#  introduced  by  for soles 

'-yTTX-  <" 
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l*t  as  utut  ptf  gt,  a,,  «|,  U|  - paraaeters  of  initial  orbit, 
and  p 4 9*  a,  e,  w - parameters  of  iatersediate  orbit  of  descent 
nodule. 

Let  us  designate  sonants  of  ties  and  polar  coordinates  of  the 
points  of  start  and  entry  of  nodule  I through  tft#  r»,  ta  and  t2»  r*, 
•i  respectively.  Obviously 

A*- f«a+A,  n<.  (2) 


•here  fu  - radius  of  planet,  daTM  - height  of  dense  layers  of  the 
ataospherau  let  us  introduce  paraaeter  p*  by  f omnia 

A-pir  (3) 

it  is  assumed  that  change  of  Ar  pglar  radius,  angular  range 


Ad  sad  duration  At  of  the  segaent  of  flight  of  nodule  I in  dense 
layers  of  the  ataosphere  to  the  point  of  finish  are  kaova  beforehand. 
Therefore  tine  t,  of  finish  and  polar  coordinates  r*,  »2  of  the  point 
of  finish  sill  be 

rt-r,— Ar, 

Oj=8a+A0. 

If  by  finish  ve  aean  the  landing  of  nodule  1 on  the  surface  of 

the  pdaaet,  then 

At* ■Am*  ^ ® 

if  the  entry  ef  this  nodule,  then 

U-hT'-h*- 0.  V-r*  »,-»►  <*-<*•  (6) 

•or  polar  radii  rt  and  rt  ve  have 
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''i-lpf+Pifi  CO*  (»1— cos  (»I (7) 
ri-P7'-\Pl+M  CO*  (*,-•)]-'.  (8) 

Characteristic  velocity  W of  initial  palse  can  be  reduced  to 
tke  fprn  [2] 


W-KAV, , 

tV^^ql+Sp'+qt—pi  — M_  —2qq  t cos  (»,  — »)— 


2S£=*>lL 


cot 


(♦i— i)}T. 


(9) 


(10) 


•here  x — Gaossian  constant,  aultiplied  by  the  square  root  of  the 
planet  mass.  For  raise  AV,  distinguished  only  by  constant  naltiplier 
from  A U,  lot  es  retain  the  nase  characteristic  velocity. 


Thrust  angle  daring  the  initial  thrust  ve  will  read  in 
reverse  notion  of  the  direction  fron  positive  transversal  to  the 
direction  of  pales*  thee  no  sill  have 


tee  - * fr »ln (»■-«)-*» «<n 

T (P\-P)  \P\+q\  co*  (»,-«,))  ' 


00 


The  signs  of  the  n a aerator  and  deaoeiaatoc  la  the  right  aide  of  (1 1) 

coincide  respectively  with  algae  sin<I>T  agd  cos  d>T- 


let  as  deter nine  angle  ® of  entry  as  the  angle  between  the 
velocity  vector  of  the  descent  aodole  and  the  plane  of  local  horizon 
at  aoaent  of  tine  ta: 

‘f* — ®e[o.-5-).  (12) 

•h*w  U,  end  — the  redial  aad  transversal  coaponents  of  velocity 
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of  nodale  I respectively  at  the  soaent  of  its  entry.  V aloes 
UnUt,Un  are  4etersine4  by  for  sal  as 

sin  (»,—).  Uk-^p-.  (13) 

*/„-*  /2^+^-p2.  (14) 

let  as  assase  at  soaent  of  tine  t2  the  orbital  nodule  has  polar 
radios  r,  and  polar  angle  t9.  Per  the  last  valaas  no  Clad 

J p <j>+q cos  v)5  + % M = J P\(P\+qi  cost/p  * ^ 

O-W+PifiCo*  (16) 

Bistanco  * botooen  nodales  at  tho  aoaeat  of  finish  is  eg  sal  to 

/-lr|+H-2r,r,cos(»,-?,)p,  (17) 

aad  co-altitude  z of  aodalo  11  at  tho  point  of  finish  sill  bo 

z—  arccos  . (18) 

2^ 


Vbe  relationships  obtained  above  show  that  as  the  basis 
parameters*  deternining  the  descent  aaneuver,  there  can  be  obtained 

P . ?.  *i.  4*  6,,  *,  (19) 

ahicb  are  dependent  aad  satisfy  the  following  connections: 

?1  cos(i >|—«)— />?—/>,?,  cos  (»,—•, )»0.  (20) 

^P'-rPq  cos  (02— «)-/>! =0.  (21) 


/>(/>+?  coif)5  J 


♦#“•1 


dv 


,-0.  (22) 


Pi(/i+fiCoi*)» 

•1— • »i-»i 

By  the  soaght  optinen  descent  naneover  an  will  naaa  the  no  soever*  to 

which  corresponds  the  least  valne  of  characteristic  velocity  AV. 


let  as  torn  to  the  calculation  of  additional  liaitations.  Let  as 
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exanine  first  only  the  liaitations  os  variables  ( esd  <l>.  hsgle  <1> 
of  entry  should  lie  in  sose  preset  interval! 

(23) 

Distaaoe  l bet seen  nodules  at  the  aoaept  of  finish  should  not  exceed 
the  fixed  aaxiaun  value  of  1: 

KL.  (24) 

Bp  introducing  additional  real  v aridities  a,  0,  conditions  (23), 
(24)  let  us  rewrite  so: 

'?*—(**••*-— '!')  (<1,~4'„„)— «J— 0,  (25) 

V,-=/-Z.+ps»0.  (26) 

thus,  with  conditions  (23),  (24)  the  problea  is  aatheaatically 
reduced  to  the  nininization  of  function  CIO)  in  the  total  of  real 
variables 

P . P.  *i.  K •.«.?.  ' (27) 

which  satisfy  conditions  (20)- (22),  (25),  (26).  Is  is  know.n,  the 
derivatives  with  respect  to  all  variables  (27)  froa  the  Lagrange 
function  should  be  equal  to  aero 

tv'+i:  **,.  <28> 

where  xt,  x2,  ....  xs  - uahaoea  constant  nu It i pliers. 

Analysis  of  the  equation  of  extrenna,  coores ponding  to  variable 
0,  shows  that  either  this  variable  is  equal  to  zero  and  condition 
i-Z.  is  fulfilled,  or  the  sought  solution  corresponds  to  the  relative 
aiaiann  of  function  AV  in  tbs  auxiliary  problea,  in  which  conditions 


(24),  (26)  are  dropped*  Therefore  it  is  possible  to  recoaaend  such  a 


* 


GE  9 

sequence  off  investigation. 

9.  The  problem  is  examined  without  taking  into  accounjt 
li ait at ions  (23) >(26).  Using  the  results  of  [ 2 ] and  the  method, 
sinilar  to  that  developed  in  article  [ 3 ]•  it  is  easy  to  show  that 
single  steady-state  solutions  in  this  case  are  the  solutions  obtained 
in  [2]  during  the  investigation  of  energetically  optiaua  single-pulse 
flight  between  elliptical  initial  and  circular  final  orbits.  The 
indicated  solutions  are  characterized  by  apsis  tangential  connections 
of  internediate  trajectories  with  initial  orbit  and  circle  of  radius 
r*.  The  trajectory,  exiting  the  apocenter  of  the  initial  orbit, 
always  reguires  less  fuel  consuaption  in  comparison  with  the 
trajectory  exiting  the  pericenter.  For  both  steady-state  solutions 
there  is  fulfilled  eguality 

«P~0.  (29)' 

In  general  case  point  (29)  lies  outside  the  interval  (23) , and  it  is 
necessary  to  change  to  the  search  for  all  relative  nininuns  of 
function  Ay,  having  kept  only  conditions  (20) -(22)  and  having  fixed 
angle  <*): 

I'-'lW  (30) 

The  analysis  and  solution  of  this  problea  with  any  fixed  angle  are 
contained  in  [3]. 

Jf  nnong  all  the  solutions  obtained  at  this  stage  there  are 


DOC  » 0674  P1GB  10 

tkosey  for  which  inequality  (24)  is  satisfied,  then  by  comparison  of 
the  characteristic  velocities  corresponding  to  then  the  sought 
solution  is  found.  In  the  opposite  case,  when  condition  (2<4)  is 
disturbed  for  all  solutions,  the  sought  optiaua  aaneuver  is 
characterized  by  equality 

/-l,  (31  > 

and  it  is  necessary  to  change  to  point  2. 

2,  A problen  is  exauined  in  which  coed it ions  (20) -(22),  (23)  are 
kept,  and  inequality  (24)  is  replaced  by  equality  (31).  For  the 
auaerical  investigation  there  can  be  used  the  aethod  developed  in  the 
following  paragraph  for  the  determination  of  the  optiaua  descent 
■aneuver  at  prescribed  values  of  / and  Q-  Interval  [Ow  <DMt]  is 
broken  down  with  the  sequence  of  points  into  a row  of  equal  parts  and 
for  each  point  there  is  calculated  the  optiaua  aaneuver.  Then  by 
comparison  of  functions  AV  there  are  fouqd  the  approximate  values  of 
angle  <D  mad  the  parameters  of  the  sought  optiaua  aaneuver.  If  the 
obtained  accuracy  is  insufficient,  there  can  be  conducted  local 
refineaent  of  these  approxiaate  values  by  one  of  the  aethods  of 
successive  approxi nations* 

he t us  assuae  now,  besides  conditions  (20) -(24),  liaitations  are 
introduced  on  variables  a and  U„.  The  jreblea  of  optiaisation  of 
function  AV.  jut  as  in  point  2,  does  net  lend  itself  to  analytical 


f| r 

f | 
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investigation.  For  its  nuaerical  solution  there  css  be  used  the 
•ethod  developed  in  the  following  paragraph  for  the  construction  of 
optima  descent  aaneuver  of  aodule  1 with  fixed  paraaeter*  l a ad  d>, 

4 which  aakes  it  possible  to  consider  the  Indicated  liaitations  on 

variables  z and  Un-  Deteraination  of  the  sought  solution  is 
perforaed  siailar  to  that  discussed  above  in  point  2*  The  distinction 
is  oahj  that  here  into  a roe  of  parts  are  broken  down  two  intervals 
[dw  Qaax].  et)—  r*.  L],  and  the  calculation  is  perforaed  for  all 

pairs  of  points,  where  one  point  is  taken  froa  the  interval  of  change 
of  angle  <t>,  and  the  other  - froa  the  interval  of  allowable  values  of 
distance  /■  The  indicated  approach  is  especially  suitable  when  high 
accuracy  is  not  reguired,  and  it  is  iaportant  to  obtain  a picture  of 
the  change  of  paraaeters  of  optiaua  aanewvers  depending  oq  the 
variations  of  vnlaaa  ‘ and  <l>. 

92*  Hatheaatical  Algorithm  of  Construction  of  Optiaua  Descent 
Baneover  Pith  Fixed  Values  of  1 nad  <l» 

Below  is  presented  a description  of  the  nethod  of  solution  of 
the  efaaiaed  problea  in  the  aasaaption  that  there  are  assigned  angle 
<i>  af  entry  of  descent  aodule  and  distance  / between  the  nodules  at 
the  tenant  of  finish,  flatheaatically  it  is  required  to  find  the 
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value  of  function  (10)  in  the  total  of  variables  (19)  * skich 
satisfy  conditions  (12),  (17) * (20)- (22),  shore  paraseters  1 «al  cl» 
are  considered  known.  The  given  aethod  aakes  it  possible  to  consider 
additionally  the  linitatiens  on  variables  z and  tin. 

Ve  will  consider  that  the  initial  orbit  is  not  circular*  where 
there  are  fulfilled  inequalities 


0<p,<Pi. 

‘ (32) 

r lu>r 2' 

(33) 

rJ»<(  + rJ» 

(34) 

rx  .>/-/*. 

(35) 

- the  distaac 

as  of  the  pericenter 

and  apocenter  respectively  cf  the  initial  orbitw  Condition  (33) 
signifies  that  the  initial  erbit  is  located  entirely  outside  the 
dense  layers  of  the  ataosphere. 

•a.  A# 

Jf  >"u<l~ri  or  ru>l+rh  then  obviously  the  descent  aaneuvers 
with  given  distance  l are  generally  ispossible.  Cases 

?i-0, 

^ii=^  ^2’  (36) 

ru-l+rt 

will  not  be  exaaiaed  here*  since  foe  each  of  thss  the  pcohles  in 

guestion  loses  its  eztrenal  character  and  is  transforned  into 
detersinate. 

let  us  change  to  discussion  of  the  aethod.  Proa  equalities  (12)* 
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(21)  >«  ob tail 

(37> 

«(*,-.)_  ''  (3>) 

where  the  signs  of  the  nnaerator  and  denoainator  in  the  last  f omnia 
coincide  with  the  signs  of  £in(%2>u)  and  cos(0s-m)  respective! y. 


Belationshlp  (13)  nakes  it  poaaihlf  to  find 


A»— 7,  *rc  cos  — ^ — , Ti-  ±1. 

o_  _ 

and  fron  equality  (20)  ee  have 


sin  (»,-*,) 


— M»fTa*:  K To! 

*V+*5 


where 


cos  (%,-«,)»  + 

Tj“±1. 

k,  — pg  cos  («•—«,)—/>! £2=/?y  sin  (u» — «»,), 

bi’=p1—p\- 


Subsequently  by  syabols  q,  lM  ltf  « we  will  aaea  the  functions 
of  pasaaeters  p,  t1(  deterained  by  relationships  (37)- (41) . lith  such 
eliaioation  of  variables  the  reaaining  unknowns  p,  1,  should  satisfy 


the  following  inequalities* 


0 <g<P, 

\n+T%-P\^2rir7, 

H+H-HX). 


Due  to  dependence  (42) . l.e»,  tie  aneaaptlea  about 

of  transfer  orbit,  ae  reduce  condition  (22)  to  the  fora 

a-0.  (45) 


eiliptieity 


where 


A— AT  H+*r  (£,-£,-#  (tin  sin  £,)J  - 

-•?  l£?)-£i’’-«i  (sin  £l’>-sin  £!«>)]. 
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I 

■here  B|,  Bs  and  E['\  Ql)-  eccentric  anonaliaa  of  aodalos  I and  II  at 
aoaeata  t1(  t,  respectively.  Tke  eccentric  anonaliaa  aca  easily 
coaputed  by  Moan  foraalaa  of  elliptic  notion.  Tina  ta-tt  of  notion 
bafoca  entry  ia  coapotad  aoi 

*2-li  - , (47) 

lot  oa  dariaa  noa  aoaa  iaegaalities  and  eatiaatioas*  abick  are 
aaad  daring  nonerical  solution.  First  of  all  fron  condition  (42)  and 

dapaadanoa 

r.  > (48) 


akara  r» 
aa  fiad 


planet ocen trie  distance  of  apocenter  of  transfer  orbit* 

P*<P^P  m 


akara 

P*  *•  - J7S-— . (50) 

/>»*=- max  j V^+PiFi  } • (51) 

Lat  oa  note  tkat  bare  the  necessary  and  sufficient  condition  (44)  of 
intersection  of  iateraediate  and  initial  orbits  is  replaced  by 
necessary  condition  (48)*  t be  ref  or  a sobseqoeatly  daring  tire  solation 
it  is  necessary  to  consider  condition  (48). 


Dependence  (43)  and  obtioos  inagoality 
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will  give 


«b«r« 


rw^  ri<  ru  • (52) 


r*<rg<r^  (53) 

r,-max  |r,„  l-rt),  r,,-min  |r„.  /+rtj,  (54) 
r»*>r»  (55) 


Th«  last  inequality  takas  piss#  toe  to  relationships  (34)f  (35). 

Condition  (53)  datarninas  two  intervals,  iq  which  lie  the  sought 

values  of  aafls  tas 


where 


•?*— •l+arcco5jt|,  J 
®f)=““j+aJ,ccosii2.  | 


7— 

« _ r» 

Mi 


IS- 


(56) . 

(57) 

(58) 

(59) 


where  obviously 

*>IS-  *<!•  l»*>— 1.  (60) 

the  developed  net hod  peraits  taking  into  calculation  the 
limitations  on  variables  z end  Um.  1st  as  assune  there  is  introduced 
requiseaeat 


(61) 


whers  *.«  - assigned  aaxiaua  value  of  co^aititade  of  aodele  XI  at 
the  ppint  of  finish  at  sonant  £.  they  fees  feraalas  (18),  (81)  ee 
obtaia 

'•>'».  (62) 

'(«S) 


f*  " K P+ri+irtl  co$zaM. 
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the  calculation  of  condition  (62)  is  reduced  now  to  conputation  of 
r.  in  dependences  (53)  -(59)  by  fornula 

r,==max  (r,.,  /— r*  /»*).  (64) 

In  a particular  case#  when  limitation  (691  is  a condition  of  direct 
visibility  between  nodules  at  the  nonent  of  finish#  distance  ru 
becones  equal  to  

r3*^  V F+n.  (65) 

Finally#  linitation  on  velocity  ee  take  in  the  fora 

Um<U,  (66) 

where  constant  0 is  assigned.  If  it  is  possible  to  disregard  the 
speed  of  rotation  of  the  at  nosphere  in  coaparison  with  value  U„, 
then  Condition  (66)  is  equivalent  to  the  xequireaent  that  the  entry 
velocity  of  nodule  I not  exceed  sone  fixed  value  of  0. 

Fron  relationships  (14)#  pi) # (#4)  ee  find 

<«d 

and#  consequently#  condition  (64)  will  fen  considered  if  inequality 
(49)  is  replaced  by 

pO<P.*- 


Subsequently#  if  we  introduce  one  or  both  liaitatioas  (62) # 

(66)  # we  will  consider  that  the  corresponding  changes  in  inequalities 
(49)#  (53)  -(59)  are  produced. 
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Let  os  foraulate  the  obtained  results.  The  probles  of 
optimization  was  reduced  to  the  search  for  paraeeter  p froa  interval 
* (49)  and  angle  A,  in  one  of  regions  (56) * (57)  so  that  conditions 

(44) * (45)  would  be  fulfilled*  and  function  AK  would  have  the 
saallest  value.  It  is  convenient  as  the  variable  during  optiaization 
to  select  paraaeter  p*  and  to  consider  corresponding  angle  as  the 
radical  of  agnation  (45)  lying  in  region  ,(56)  * (57)  * for  which 
condition  (44)  is  fulfilled. 

Depending  on  the  concrete  selection  of  paraaeters  yx,  ya  four 
types  of  aaneovacs  should  bo  investigated,  hot  as  introduce  value 

j-  5— ~ . 7i=~  ±1.  T*“±l.  (69) 

which  takes  values  j ■ 1*  2#  J#  4 for  thuse  types,  she  descent 
trajectories  for  j = 1*  3 (y«s1)  differ  froa  the  tra jectories  for  j * 
2*  * *1)  by  concrete  selection  of  the  point  of  start  in  one  of 

the  tao  points  of  intersection  of  initial  and  interaediate  orbits. 
Therefore  in  one  case  the  flight  path  of  aodule  I will  eaerge  on  the 
initial  aegaeut  froa  the  region  liaited  by  the  initial  orbit*  and  in 
the  other  case  will  entirely  lie  inside  the  indicated  region. 

Further*  aaneuvers  j ■ 1*  2 (lt*t)  differ  froa  aaneavers  J » 3*  4 
(!*»  *1)  by  the  fact  that  at  the  aoaent  of  finish  in  the  first  case 


the  orbital  aodule  leads  the  descent*  i.e«*  has  large  polar  angle.  In 
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tli«  second  c«N(  conversely,  the  polar  angle  of  the  descent  sod  ale  is 
larges  than  t|e  polar  eagle  of  the  orbital.  Let  ns  note  tlrnt  sith 
each  eoa  peri  see  angles  di.  A,  shonld  be  brought  to  interval 

I «!-«»! 

Bet  ns  present  a brief  description  of  a coapoter  prog  ran,  using 
the  algoritha  indicated  above. 

I.  A search  is  conducted  for  the  approx in ate  desired  value  of 
paraaeter  p bf  aeaas  of  the  total  survey  of  interval  (49)  r vhich  is 
broken  down  into  a egual  parts.  For  each  point  of  division  there  are 
found  all  the  radicals  of  agnation  (45) , vhich  lie  in  regions  (56) , 
(57)  and  for  vhich  is  fulfilled  condition  (44) • During  the  coaputiag 
of  radicals  there  are  taken  rough  constants  of  accuracy,  vhich 
provides  the  quickness  of  operation  of  this  prograa  unit.  Proa  the 
aultitude  of  all  the  sorted  out  values  of  paraaeter  p and  the 
cadicbls  corresponding  to  then  there  is  selected  pair  p,  h)t  vhich 
corresponds  to  thn  saallest  value  of  function  AV.  These  quantities 
are  taken  as  appro xisate  options  values. 

II.  There  is  found  the  exact  opt in on  value  of  variable  p and 
other  paraaeters  of  the  sought  naaeuver*  Let  ns  note  that  quantity  AT 
ns  a function  of  paraaeter  p vith  opt in on  selection  of  radical  t,  can 
have  discoatikuities  of  the  first  type  and  regions  of  uaiaportaace. 
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■here  090a t ion  (45)  does  not  have  radicals  of  the  required  type  at 
all.  Bae  to  the  indicated  special  character  of  the  problem  it  is 
expedient  to  refine  the  optima  values  of  pare  asters  until  the 
obtaining  of  certain  accuracy  by  the  aethod  of  successive 
ap pros iaat ions.  In  this  case  in  each  approxiaation  the  interval 
between  two  values  of  paraaeter  pf  adjacent  to  the  options  value  of 
the  preceding  approxiaation*  is  broken  down  intc  sons  preset  nunber 
of  pasts  and  calculation  is  perforaed  for  all  the  points  of  division. 
By  coe pari son  of  the  corresponding  characteristic  velocities  the 
optiaaa  value  of  paraneters  of  the  given  approxiaation  is  deterained. 
During  coa potation  of  the  radicals  of  function  a there  are  taken 
constants#  providing  the  assigned  accuracy  of  ccaputations. 

In  conclusion  let  us  elaborate  on  the  procedure  of  conputation 
of  radicals  of  function  (46)  for  any  fixed  value  of  paraaeter  p.  For 
deterainacy  we  are  Halted  by  descent  aaneuvers*  for  which  the 
overall  tine  tr-t,  of  the  ease  over  is  exactly  less  than  thre  period  of 
aotioa  along  tto  Initial  orbit.  In  accordance  with  this  we  will 
coapate  the  eccentric  anoaalies*  entering  the  equalities  (.45)- (47)* 
taking  into  accooat  conditions 

Ei 

£?>+2*>£y>+2*-*>£y>>£w. 

where  6 - any  assigned  saell  gaeatlty. 


legions  (56)*  (57)  are  divided  into  I eg sal  parts  and  f auction  A 


DOC  * 0674 


PAGE  20 


is  coaputad  successively  for  all  the  points  of  division. 

Let  us  indicate  one  special  feature  of  the  given  nethod.  For  any 
value  of  angle  4,  first  there  is  deterniaed  the  geonetric  picture  of 
the  aaneuver,  i.e.  , nutual  location  of  orbits,  and  then  there  are 
selected  the  specific  legs  of  flight  of  the  nodales  along  orbits  fron 
the  cpndition  of  fulfillnent  of  inequalities  |?Q).  Therefore  critical 
values  of  angle  A,  can  exist,  which  correspond  to  discontinuity  of 
function  St-E|«  if  switching  of  the  leg  of  flight  occurred,  on  the 
descent  orbit,  and  functions  B^-Bp,  if  switching  occurred  on  the 
trajectory  of  notion  of  nodale  XX.  At  the  points  of  intersection  one 
of  the  liniting  values  of  the  corresponding  functions  is  equal  to 
zero,  and  the  anount  of  discontinuity  is  equal  to  2v.  It  is  easy  to 
see  that  the  radicals  of  function  A can  lie  only  at  finite  distance 
fron  the  critical  points.  Let  as  select  aeaber  I so  that  the 
variation  of  functions  Ex—E,.  £y>—  £{'»  op  each  segaent  of  subdivision 
of  intervals  (56),  (57)  foes  not  eaceed  seee  constant  A.  The 
variation  of  function  on  the  segaents  containing  critical  points  will 
then  he  not  less  than  2*-A.  By  increasing  I and  selecting  constant  A, 
it  is  always  possible  to  achieve  falfilleent  of  condition  A<2r-A  and 
thereby  obtain  criterion,  which  nates  it  possible  by  the  eagnitude  of 
fluctuation  of  function  A to  judge  whether  switching  was  inside  the 
considered  segaent  or  not. 
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let  us  return  to  the  description  of  the  procednre  of  computation 
of  radicals.  If  on  the  boundaries  of  son#  segneat  of  division  of  the 
intervals  (56),  (57)  function  A has  different  signs,  and  its 
flucteation  is  less  than  constant  i,  then  within  this  segneat  is 

« 

found  the  radical  of  function  A , which  is  coaputed  with  the  aid  of  a 
series  of  interpolations. 

• 

let  us  elaborate  aore  os  one  specific  feature  of  the  algoritha. 
If  for  sose  value  of  condition  (44)  is  disturbed,  then  we  will  say 
that  this  value  is  found  in  the  region  of  nonexistence  of  function  A. 
■hen  during  sorting  out  of  points  fron  intervals  (56),  (57)  we  fall 
in  the  region  of  nonexistance  of  function  A,  then  b;  successive 
division  of  the  given  segaent  in  half  we  find  the  boundary,  of  the 
region  of  nonexistence.  Then  we  investigate  the  region,  one  of  the 
boundaries  of  which  is  the  last  point  during  sorting  out  of  points  of 
subdivision  of  intervals  (56),  (57),  and  the  other  - the  boundary  of 
the  region  of  nonexistence,  aith  the  fulfilment  of  conditions  of  the 
presence  of  radical  in  the  indicated  region  the  radical  is  cdsputed 
by  a series  of  interpolations.  An  analogous  procedure  is  pertforsed 
when  fros  the  region  of  nonexistence  functions  A transfer  to  the 
values  of  angle  ts,  for  which  condition  (A 4)  is  fulfilled. 

After  all  the  radicals  of  function  A are  found  for  the  paraaeter 
p is  question,  fros  then  is  selected  a radical,  which  the  snallest 
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value  of  function  AT  corresponds  to.  iith  this  the  deteraination  of 
the  radicals  of  function  A is  finished. 

i * 

§3.  laaerical  Bran pie 


As  an  eraaple  of  the  use  of  the  algor itha  developed  in  §2  let  os 
exaaiae  the  prohlea  of  sending  fron  a space  vehicle,  aovinjg  along 
elliptical  orbit  of  a Tenus  satellite,  probe  (nodule  I)  for 
investigation  of  the  upper  layers  of  the  Tenus  ataosphere.  Re  will 
ainiaize  the  fuel  consunpticn. 


ts  of  the  initial  orbit  of  the 


space  vehicle 


10000  km,  *.=0.28, 


The  hhight  of  periceater  of  the  orbit  abose  the  surface  of  Tenus  vill 


be  1000  ha,  and  the  height  of  apocenter  - 6600  ka 


of  constant  I and  r#  lot  no 

**~33423  10»  KjitjcM*, 
r2  =6200  KM. 

the  toooat  of  finish  coinci 


aoaent  of  entry,  t 
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Bet  as  request  that  at  the  soseat  of  finish  values  l mmk  had 
assigned  values  and  the  condition  of  direct  visibility  bee seen 
nodules  vas  fulfilled- 


Calculations  sere  conducted  on  an  1-20  coaputer.  Two  values  of 

distance  were  selected 

/« 1200  KM,  /— 1600  km  (74) 

and  the  interval  af  ctaefe  ef  an  fie  • was  ezaaiaed  with  pulse  duty 

(actor  5*  froa  aero  to  30°-  The  results  of  calculations  are  presented 

la  Tables  1 aad 


b) 


- 

l *J 

— rrr 

W.  km  fee* 

<1 — 7|, 

/ ft* 

t Cl 

# 

/ 

K,  KM 

« 

•i 

•* 

*» 

at 

*t 

X 

(/.*,  KM/cetc 

O’ 

3 

9497,4 

044V? 

194*99 

368, *61 

2.95 

8*61 

0,348 

126*14 

5305 

36*03 

8,394 

5 

3 

9170.7 

04343 

222.83 

355.63 

350.25 

15,74 

0,329 

149,09 

3259 

.14,32 

8,321 

10 

3 

6949,2 

ojun 

22842 

339,73 

316.40 

19,59 

0,401 

156,54 

2787 

20,80 

8,268 

15 

1 

8394,2 

04616 

24544 

330,91 

311,27 

31,62 

0,601 

163,91 

2055 

2,21 

8,122 

20 

2 

7486,4 

047*2 

2*1.72 

33142 

33160 

56,05 

1480 

195,98 

1045 

14,73 

7,828 

» 

> 

39141 

mj o 

31548 

73.93 

1,773 

214,19 

645 

19,83 

7,634 

» 

t 

«7UM 

MW 

9114* 

334,02 

33740 

85,84 

2,480 

1 

222,28 

468 

22,96 

7,526 

• » 1200  Its. 


■sys  Ms)  ks-  |b)  ka/s-  <c>  a 
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D 

D 

rngLi 

aa 

a 

R9 

a 

a 

0 

a 

U — it* 
cett 

B 

0° 

3 

9298.9 

0,3333 

219,64 

372.46 

1 1 

361  >4 

12,46 

0,268 

156.68 

3589 

56,84 

8,350 

5 

3 

9321,2 

0.3448 

214.01 

350,65 

340,66 

10,29 

q.261 

167,25 

3491 

52,26 

. 8,.’55 

10 

3 

9172.4 

0,3633 

212,90 

333.36 

325.39 

11,92 

0,328 

167,40 

3276 

40,44 

8,321 

15 

3 

8888.3 

0.3903 

216.58 

319.84 

315,04 

16,38 

0,453 

106.52 

2918 

23,71 

8,253 

20 

1 

8459,2 

0.4242 

222.85 

310,42 

310,42 

24.15 

0,655 

102,30 

2553 

1,02 

8,140 

25 

2 

7696.6 

0,4560 

249.41 

300.20 

312,42 

41.54 

1.003 

180,46 

1580 

15,80 

7,904 

30 

2 

702W 

0,5103 

270,74 

310,54 

3ISJ9 

50,02 

1,552 

197,20 

999 

24,99 

7,645 

HU«  2.  rtructora  of  optim  •««•«««  foe  1 ■ 1400  ka. 


Key:  sua  as  Table  1. 


Let  os  note  that  foe  both  *ala«s  of  distance  l the  aiaiaoa  of 
characteristic  velocity  as  a foactioa  of  angle  4>  is  reached  for 
soaa  v aloes  of  angle  ^€(0*^),  shore  this  velocity  rapidly  rises  for 
d>>10+*  It  is  interesting  to  note  also  that  velocity  for  angles 

0>,  net  exceeding  10-15°  for  optinan  aaneuvers,  barely  depends  on 
asgle  ® and  is  egaal  to  approxinatel;  8*3  kn/s. 


i 
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